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A number of evidences suggests that thick-film resistors are close to a metal-insulator transition and 
that tunneling processes between metallic grains are the main source of resistance. We consider as 
a minimal model for description of transport properties in thick-film resistors a percolative resistor 
network, with conducting elements governed by tunneling. For both oriented and randomly oriented 
networks, we show that the piezoresistive response to an applied strain is model dependent when 
the system is far away from the percolation thresold, while in the critical region it acquires universal 
properties. In particular close to the metal-insulator transition, the piezoresistive anisotropy show 
a power law behavior. Within this region, there exists a simple and universal relation between the 
conductance and the piezoresistive anisotropy, which could be experimentally tested by common 
cantilever bar measurements of thick-film resistors. 
PACS numbers; 72.20.Fr, 72.60. -hg, 72.80.Ng 



I. INTRODUCTION 

Thick-film resistors (TFRs) are a particular class of 
granular metals made of Ru02, Bi2Ru207, or other 
metal-oxide granules embedded in an insulating glassy 
matrix. The importance of TFRs in applicative prob- 
lems is due to their large piezoresistive response which, 
together with high stability and low manufacturing costs, 
makes the. TFRs highly competitive as pressure and force 
sensors. H'El Hence, the full understanding of the factors 
governing the piezoresistive response is a very important 
issue in optimizing the performances of TFRs. The inter- 
est in the piezoresistive response however concerns also 
more fundamental aspects connected, as we discuss in 
the following, with the role of transport anisotropy in 
disordered systems. 

In practical measurements, the piezoresistive response 
is obtained by the so-called gauge factor commonly 
defined as K = SR/eR, where SR is the change of resis- 
tance R for a given applied strain e. When the strain 
is directed along the voltage drop, the change in resis- 
tance determines the longitudinal gauge factor Kj^. Usu- 
ally, in TFRs the gauge factor ranges from Ki^ 2 up 
to K]^ 35 or more,0 following roughly the empirical 
relatior 
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Both the high values of Ki^ and Eq. (|^) point towards a 
scenario in which tunneling processes behween metallic 
grains are the main source of resistance.QO In fact, if d 
is the mean tunneling distance, then transport between 
adjacent metallic grains is approximatively 



R oc exp(2d/0. 



(2) 



where f is the localization length. Let us assume that 
under an applied strain e the distance d changes to d ^ 



d + Sd = d{l + e), where we have set e = dd/d. The 
strain- induced change of R is therefore SRJ.R ~ 2de/^, 
which implies the empirical relation Eq.(|l|).H 

In obtaining this simple result, we have assumed that 
the direction of intergrain tunneling is along the direc- 
tion of the applied strain. In this way we obtained an 
upper limit of the longitudinal gauge factor Ki^. In struc- 
turally isotropic systems, however, the direction between 
two neighbouring grains is oriented randomly and can 
be even ortogonal to the direction of e. This leads to 
a reduction of Kj^ and to a contemporary appearance 
of a nonzero transverse gauge factor Kt defined by the 
change of resistance when the voltage drop is ortogonal 
to the direction of e. In addition, if we consider also that 
TFRs are disordered conductors, then it is natural to ex- 
pect that, in fiowing from one end to the other of the 
sample, the microscopic currents are subjected to dras- 
tic changes of direction sampling therefore regions where 
the intergrain distance is d{l+e) and regions where d re- 
mains unstrained. This effect contributes to the decrease 
of Ki^ and to the enhancement of K^B 

What is commonWppbserved in TFRs is that i^L and 
Kt are comparableJ3l3 indicating therefore that the cur- 
rent tortuosity induced by disorder and random bond 
orientations is an important element and needs further 
investigation. To this end, we have recently addressed 
the problem of the piezoresistance in a simplified model 
of TFRs (or more generally of granular metals) consist- 
ing of a random resistor network with a fraction p of 
bonds having resistances of the form of Eq.(|^) ancLthe 
remaining fraction I — p having zero conductance.Ei In 
accord with the qualitative discussion of above, we have 
found that as p is reduced from the unity, the longitu- 
dinal gauge factor is diminished and the transverse one 
is enhanced. This trend holds true all the way down to 
the critical bond concentration pc at which the conduc^ 
tance G vanishes and Ki^ — s- Kt for an uniaxial strain.El 
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The observation that Ki^ ^ Kt in TFRs could therefore 
be an indication that these systems are close to a metal- 
insulator transition, an hypothesis sustained by the very 
low metallic volume fractions v in TFRs and the observa- 
tion that the conductance G follows a power law behavior 
G ~ (u — Wc)*, where Vc jf Ae critical volume fraction and 
t is a critical exponent .Lril2l 

The aim of this paper is to study further the prob- 
lem of piezoresistive response in random resistor net- 
works by focusing on the equivalence between transport 
in anisotropic networks and the piezoresistive response 
of isotropic ones. Close to the percolation thresold, the 
two systems are equivalent in the sense that transport 
is governed by the same critical exponents. In addition 
we consider also the effect of random bond orientations 
and provide explicit formulas from which it is possible to 
extract the piezoresistive anisotropy. 



II. PIEZORESISTIVE ANISOTROPY OF 
RANDOM RESISTOR NETWORKS 

In this section, we describe in full generality how to ex- 
tract from a random resistor network model the relevant 
piezoresistive coefficients. Although real materials like 
TFRs are quite complex from the point of view of both 
microscopic electronic processes and microstructure, we 
believe that the essential physics of the piezoresistive re- 
sponse is already contained in our simplified model. Pos- 
sible limitations of our approach are however discussed 
in Scclvl. 



In common piezoresistive measurements, the sample 
is subjected to a geometrical distortion characterized 
by a strain field with coefficients en {i — x, y, z) and 
the piezoresistive response is obtained by recording the 
change of resistance (conductance) under the effect of 
£ii. In the following, we shall consider only situations 
in which the shearing strains can be neglected, so that 
Eij — Q ior i ^ j. In the absence of imposed strains, we 
assume that the sample is isotropic and characterized by 
a conductance G. Hence, under the effect of e^, the con- 
ductance for a voltage drop in the i direction changes to 
Gi = G -\- SGi, where the variation SGi can be expressed 
in terroa of the conductivity change Sai and a geometric 
factor :Ej 



SG^ 
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where a is the conductivity for the unstrained sample 
and the indexes j, j, and k assume the values x, y, and 
z with cyclic permutations. 

We are interested on the intrinsic conductivity change 
Sai which is governed by the microscopic electronic pro- 
cesses taking place in the bulk. In agreement with the 
empirical relation Eq.(^), we assume that bulk trans- 
port is mainly governed by tunneling processes between 
metallic grains and we consider a cubic random-resistor 



network whose unstrained bond conductances are either 
proportional to exp(2(i/^), where d is the distance be- 
tween two neighbouring sites, or zero. The probability 
of having zero bond conductance can be independent of 
the particular bond considered, as in bond percolation 
models, or correlated, as in site percolation models. The 
following considerations apply also for more complicated 
distributions, like the ones arising in segregated site per- 
colation models where all the sites within spheres of a 
given radius centered at random are excluded from the 
network, E2I or when the non-zero bond conckiGtances have 
some distribution of tunneling probability.H'lllI 

To model the effect of the applied strain, let us assume 
that the resistor network is embedded in a homogeneous 
elastic medium and that the elastic coefficients of the 
network and the medium are equal. Moreover we assume 
also that the directions of the bonds in the network are 
aligned to the x, y, and z axes. Therefore, in general, 
when a strain field Sa is applied on the sample, the con- 
ducting bonds change their values according to their ori- 
entation with respect to en. For example, the tunneling 
distance for a bond directed along the i axis changes to 
d di = d(l -|- en) which, up to the first order in en, 
implies a bond conductance change 
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Now if we consider an uniaxial strain along, for example, 
the X direction {exx = £, £yy = £zz — 0), thus the bond 
conductances directed along i = x are modified according 
to Eq.(Q), while those directed along the y and z axes 
remain unchanged. Therefore, up to first order in e, the 
resulting conductivities can be expressed as ax — <J — 
crFije, and ay = a^ = a — aT±e, where we have defined 



Fii = 

r, = 
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as the longitudinal |-ajid transverse piezoresistive coeffi- 
cients, respectivelyJld The above reasoning holds true 
also for uniaxial strains along the y and the z axis and, 
since the problem is linear, for a general strain field en 
(i — X, y, z) the conductivity variations Sai reduce to: 
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As already pointed out in the introduction, the com- 
monly measured quantities used to extract the strain sen- 
sitivity of the sample are the piezoresistive gauge factors 
which for a general strain are defined as: 
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If we assume that the strains Su are known, then 
Eqs.(||j7|) permit to express the different gauge factors 
Kij in terms only of the two intrinsic piezoresistive coef- 
ficients r|[ and r^. For example, in a typical cantilever 
bar experiment with the cantilever main axis directed 
along the x direction, the strains are approximatively 



e, Eyy — —ve, and £z 



-v'e. where v and 



are the Poisson ratios-|p£| the cantilever and the resistive 
sample, respectively .IIjO By using Eqs.(H.M), the longi- 



tudinal and transverse (Kt) piezoresistive gauge 

factors are: 

ifL = i^xx = (i + r,|) + (i-ri)(z/ + z/'), (9) 
ifx = = -(1 + r||)i/ - (1 - ri)(i - v'). (10) 

From a measurement of i^L and -f^T, the above expres- 
sions permit to extract Fn and F^ which are intrinsic 
quantities, i. e. they do not depend on the particular 
strain applied. Note that instead of measuring Ki^ and 
i^T, the piezoresistive coefficients can also be extracted 
from Eqs.(m^) by measuring a gauge factor under two 
different imposed strains. For example, as in Ref. 
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cantilever bar measurement provides a first value of A'l, 
Eq. , and a successive measurement of the longitudinal 
gauge factor under hydrostatic pressure provides a sec- 
ond value of ■ These two values are then sufficient to 
extract F|| and F^. 

The importance of F|| and V± resides on the fact that 
they permit to extract useful informations on the per- 
colative nature of transport. To illustrate this point, we 
show in Fig. 1 the results of numerical Monte Carlo cal- 
culations of F|| and F^ (filled circles) for a site perco- 
lation model in which a concentration of sites x is re- 
moved at random. For each missing site the bond con- 
ductances connecting the six neighbouring sites are set 
equal to zero, while the remaining bonds have conduc- 
tance exp(2(i/^) with 2(i/^ = 4 when unstrained. Note 
that in this model, the fraction of conducting bonds p is 
equal to x^, since a bpad is present only if the sites at 
both ends are present .lij We calculate the total conduc- 
tance by solving numerically the Kirchhoff equations for 
all the nodes of the cubic network. In practice, we im- 
pose a unit voltage difference between two opposite sides 
of the network with periodic boundary conditions to the 
remaining sides. The piezoresistive coefficient F|| (Fj^) is 
then obtained by calculating, for a fixed configurations of 
bond resistors, the difference in conductance when e = 
and e = 0.01 for strain directed parallel (ortogonal) to 
the imposed voltage drop direction. 

When there are no missing sites (x = 1, p = 1), the 
current flows exclusively along paths directed along the 
direction of the voltage drop. In this case, the longitudi- 
nal piezoresistive coefiicient is F|| = 2c?/^ — 4 while the 
transverse one Fj^ is zero [see Eqs.(|,|)]. When sites are 
removed {x < 1), Fy gets reduced and at the same time 
the transverse coefficient is enhanced in such a way that 
Fy > F_L. This behavior is due to the fact that as sites are 
removed from the network, the missing bonds force the 
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FIG. 1. Monte Carlo calculations of longitudinal (Fy) and 
transverse (Fx) piezoresistive coefficients of site percolation 
cubic resistor networks with number of sites up to 40 x 40 x 40. 
p = a;^ is the fraction of bonds present for a concentration x 
of sites. The critical bond fraction is pc — ~ 0.098. Filled 
circles refer to networks with bonds oriented along the x, y, 
and z axes, while the open circles are the effective longitu- 
dinal (Fii) and transverse (Fx) coefficients resulting from an 
average over bond orientations [see Eqs.( p^ , p3| )]. 



current to flow also along directions perpendicular to the 
voltage drop. Hence, when the strain is along the volt- 
age drop, the current visits also regions where instead the 
bonds are unstrained (reduction of F||) while when the 
voltage drop is perpendicular to the direction of e, the 
current is influenced also by regions where the bonds are 
strained (enhancement of F^). This trend gets amplified 
as X is further reduced and at the percolation thresold 
{pc = xl c± 0.098) F_L — > F||. Quafitatively, this behavior 
is observed also when the bond conductances have some 
distribution or when different statistics like that of bond 
percolation models are considered.alliJ Hence, as a gen- 
eral rule, the longitudinal and transverse piezoresistive 
coefficients become equal as the random resistor network 
reaches its percolation thresold. 

In relation to the gauge factors A'l and Kt, it is worth 
to point out that, from Eqs.(|9|,p^, JsTl — Kt = (F|| — 
Fj^ + 2)(1 + i/) in accord with Ref. However, since 
F|| > F^ for Pc < H-^ 1, it is not justified to assume 
i^L — Kt — 2{l + iy),i2 unless the system is very close to 
its percolation thresold. Instead a more general relation 
is Ki^ — Kt > 2(1 -I- i^), which holds true in the whole 
range olp>Pc values. 

Summarizing the results of Fig. 1 , in the presence of an 
external uniaxial strain, the network is anisotropic when 
p \ while is electrically quasi- isotropic when p ^ pc- 
The metal-to- insulator transition can therefore be viewed 
also as a piezoresistive anisotropic-to-isotropic transition 
when p Pc- We call the quantity which characterizes 
this transition the piezoresistive anisotropy factor x, de- 
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fined as: 



rii 



(11) 



As shown in Fig. 2, x (filled circles) is equal to the 
unity for p = 1 while it decreases monotonically as the 
system moves towards its percolation thresold and, as 
shown in log-log plot of the inset of Fig. 2, x goes to zero 
by following a power law behavior in the critical region 
{p~Pc) < 1: 

X-{p-Pcf. (12) 

where A is a critical exponent A fit of the numerical 
data of Fig.2 to Eq.(|l|) leads to A = 0.44±0.07. We have 
tested that Eq.(|lJ) holds true also for bond percolation 
models and for conductiiig bonds having a distribution 
of tunneling probability^ 



The power law behavior of Eq.(12) is directly related 



to the transport properties of anisatxopic bond percola- 
tion models studied some time ago.E3lla In those works, 
percolating networks with random bond conductances 
were defined in such a way that the conducting bonds 
were equal along two directions, for example y and z, 
but different from those along the third direction, that is 
X. Topological considerations renownalization group 
analysis and numerical calculationala showed that the 
quantity Oyjax — 1 — oiip) has a power law behavior 
Q^(p) ^ {p ~ Pc)^ close to the percolation threspM, inde- 
pendently of the microscopic bond anisotropyllZI In our 
model, bond anisotropy is induced by the effect of an ap- 
plied strain, Eq.(^, so that the two situations are equiv- 
alent. It is in fact easy to see from Eqs.(^,p|) that for our 
piezoresistive model ay /ax — 1 ~ (Fn — F_l)£- oc x^ so that 
the two exponents A and A' are equal. Since the exponent 
A does not depend on the anisotropy of the microscopic 
bond conductances, a measurement of the piezoresistive 
anisotropy factor x is in principle a practical and alter- 
native way to study the role of anisotropy in percolative 
conductors. 

In relation to TFRs, a measurement of the piezore- 
sistive anisotropy factor % as a function of metallic vol- 
ume concentration would be an useful tool to character- 
ize their transport properties. The advantage of x over 
other quantities like for example Ki^ — K^^ is that it mea- 
sures the anisotropy without being seriously influenced 
by other factors like the change of tunneling distances 
with the metallic concentration x. By lowering x in fact 
the mean distance between metallic granules should in- 
crease leading to an enhancement of the tunneling dis- 
tance d. This situation can be modeled by assuming that 
d is some function of the bond probability p, d{p) — df{p), 
where f{p) increases as p is reduced and is independent 
of strain. It is easy to see that, if the unstrained conduct- 
ing bonds have conductance exp(2d/^) = exp(2d/(p)/^), 
then the longitudinal and transverse piezoresistive coeffi- 
cients get multiplied by f{p)- This factor exactly cancels 
out in the definition of Xj Eq. ([ll|) . 
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FIG. 2. Piezoresistive anisotropy factor as a function of 
p ~ Pc for the data of Fig.l. Filled circles (x) refer to the x, 
y, z oriented network while open circles (x) are the results of 
an average over bond orientations [see Eq.(^)]. In the inset, 
the same data are plotted in a log-log scale. 



Let us conclude this section by pointing out a relation 
between conductance and piezoresisitve response differ- 
ent from Eq.(^. As already pointed out in the intro- 
duction, TFRs should be quite close to a metal-insulator 
transition, since the volume concentration is small [js4 
the conductance G shows a power law behavior .LTtHl 
Hence, if TFRs are actually in the critical regime, then 
as a function of metallic concentration the piezoresistive 
anisotropy factor should follow Eq.(|l^) and consequently: 
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where t is the critical exponent of the unstrained conduc- 
tance G ^ {p — PcY- 



III. AVERAGE OVER BOND ORIENTATIONS 

Equation (|^) of the last section has been obtained by 
assuming that the bond directions are aligned parallel to 
the X, y, and z axes of an ortogonal frame. In this way, 
under the effect of an applied strain en, the conducting 
bonds change according to the simple relation of Eq. (||) . 
In this section we generalize Eq.(|7|) in order to describe 
networks with bonds oriented at random with respect to 
the reference frame. As already pointed out in the intro- 
duction, we expect that in this more realistic situation 
the longitudinal and transverse piezoresistive responses 
should be different from those obtained for the oriented 
network, at least for p ~ 1. 

Instead of studying a network with bond directions 
changing from site to site, we consider an ensemble of per- 
fectly cubic networks oriented at random. The piezore- 
sistive responses will then be obtained by averaging over 
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all network orientations. We start by choosing the three 
bond directions of the cubic network as oriented along the 
x', 2/', and z' axes which constitute an orthogonal frame 
A' rotated by R with respect to the reference frame A 
specified by the x, y, and z axes. The rotation matrix R 
is given by the triple product of successive rotations de- 
fined by Euler angles. To fix an explicit rotation matrix, 
let us consider a z-x-z transformation given by a rotation 
of angle (p around the z-axis, followed by rotations of 9 
and ijj around the the new x and z axes, respectively. 
With this choice of angles, ^an take values from to tt, 
while (f> and tp from to 27r.ll3 By using the usual rules of 
strain transfornmtions, the strains e'^ in reference frame 



A' are given 



by:^ 



^ij — ^ RiiRji^u, 



(14) 



where Rij are elements of the matrix R and, as before, 
we have neglected the shearing strains with respect to 
frame A. Now, let us consider in the absence of strains 
a conducting bond directed along the x' direction having 
tunneling length dx' = d. Under the effect of Eg. (14), 
the bond length changes to 

^{d + Sdx')^ + {Sdy^)^ + iSd,,)^ 

= dy'(l+4j2 + (,,^)2 + (,,J2 

= d(l + e;j + 0(4J + 0(4,) + 0(4J, (15) 

where Sdi' is the variation of length along the i' = 
x',y',z' axis. Compared to e'^x, the effect of the shear- 
ing strains e'^y and e^^ is therefore of higher order and 
can be neglected. Hence for a bond along the i' di- 
rection the bond tunneling length changes simply to 
d ^ dii — d{\ + e--). At this point, with respect to 
the A' frame, we can follow the same reasonings of last 
section so that the conductivity changes are of the same 
form of Eq.(0): 



= -rj.A-Kr^-r||)4, 



(16) 



where 5a[ is the variation of conductivity a[ in the A' 
frame and we have introduced for later convenience the 
volume dilatation A — Y^- e'jj . 

Now, the local Ohm law is j- = cj[E[, where j[ and 
E[ are the components with respect to A' of the cur- 
rent density and the electric field, respectively. These 
are related to the corresponding quantities in frame A 



ji=T,iRijj3 andi;,' 



where: 



RijEj, so that ji = CT. 



(17) 



Since the dilatation A is invariant under rotation, equa- 
tion drn) reduces to: 



i^ij 



%(7(i-riA) + (7(rj 



(18) 



where we have used Eq.(|l^) and 6ij is the Kronecker 
delta symbol. Finally, the average over the ensemble of 
networks with different bond orientations is achieved by 
performing the average over the angles 9, </>, and ip ap- 
pearing in R. For a general function f{9,(j),'ip), this is 
achieved by 



/; d9 /o'" dcP J^^ d^J{R)f{9, </), ^) 
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d(f> dil;smi9)f{9,^,ib), (19) 



where J{R) = sin(6')/8 is the Jacobian of R. It is easy 



to see that when e'n in Eq.(|l8[) are replaced by Eq.(|l4| 
all the off-diagonal components of vanish under an- 
gle integrations, (cry ) = Sij{(7i), and that the averaged 
conductance variations {Sai) = (ui) — cr reduce to: 



i^ = -r,A + (r,-r„)^ 
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By using A = Sjj , this expression can be recast in 
form more similar to Eq.(0): 



r II Eli - rj_(£jj + £fcfe), 



(21) 



where we have introduced the effective longitudinal and 
transverse piezoresistive coefficients 



SFii +2rx 



r, = 



4ri 



(22) 
(23) 



If all bonds are present and have all equal conductances 
exp(2(i/^), then from the above equations we obtain 
f|| = 3r||/5andf_L = Fn/S, since r_L = OandFn = 2d/i. 
Therefore, as expected, with respect to a network with 
bonds oriented along the x, y, and z axes, the longitudi- 
nal response is diminished and the transverse one is en- 
hanced. On the other hand when the network approaches 
to its percolation thresold, Fn F^, also the effective 

responses become equal F|| — > Fj^. This behavior is sum- 
marized in Fig. 1 where the open circles refer to F|| and 

F_L obtained from the corresponding coefficients for the 
oriented network via Eqs.(22 2^). 

The effect of averaging over all network orientations 
is therefore substantial away from the critical region and 
unimportant whenp— pc ^ 1. This reflects the universal- 
ity of the piezoresistive response close the metal-insulator 
transition which is more evident in the behavior of the 
effective piezoresistive anisotropy factor 
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r,i - r. 



2X 



3 + 2(l-x)' 



(24) 



which, according to Eq.(n2|), shows the power law behav- 
ior X — 2x/5 ^ {p — Pc)^o'i: p — Pc ^ 1 (open circles of 
Fig. 2). 

IV. DISCUSSION AND CONCLUSIONS 

In the previous sections, we have shown how to ex- 
tract the intrinsic piezoresistive responses Fn and r± of 
a granular metal from the knowledge of the gauge fac- 
tors and the values of applied strains. The so-obtained 
longitudinal and transverse coefficients permit to define 
a piezoresistive anisotropy factor x. Eg. ([TTI) , which fol- 
lows a power law behavior if the system is sufficiently 
close to its percolation thresold. Concerning the prop- 
erties of TFRs, a measurement of x as a function of 
the volume fraction could be a useful tool to investi- 
gate the closeness to a metal-insulator transition. We 
have also considered the effect of random bond orienta- 
tions by employing an average procedure over an ensem- 
ble of networks with different orientation. This scheme 
has permitted us to obtain a straightforward and simple 
generalization of the results obtained for networks with 
fixed orientations. Unfortunately, we are not aware of 
published data from which it is possible to extract the 
piezoresistive coefficients as a function of metallic con- 
centrations, so our predictions cannot yet be tested. We 
hope that the present study will encourage acquisition 
of data capable of estimating x as a function of metallic 
concentrations. 

Before concluding, it is useful to remind possible limi- 
tations of the theory presented in this work. First of all, 
all through the previous analysis we have assumed that 
the granular metal is elastically homogeneous, so that the 
strain values inside the bulk can be considered as inde- 
pendent of the position. However, TFRs are strongly 
heterogeneous composites in which quite stiff metallic 
granules are embedded in an elastically soft glassy ma- 
trix. This gives rise to important fluctuations in the lo- 
cal values of the strain. Since the microscopic tunneling 
processes are affected by the local rather than the macro- 
scopical strain values, such an elastic heterogeneity can 
influence the piezoresistive response in an important way. 
However, this effect influences the absolute values of Fn 
and F^, but should leave the piezoresistive anisotropy 
factor X relatively unaffected. This is indeed confirmed 
by the analysis Cjf_.a simplified model of elastic hetero- 
geneity in TFRs, til from which actually it results that 
(Fy — F_l)/F|| is independent of the metal and glass elas- 
tic coefficients. This is of course only a preliminary result 
which needs further investigations. 

Another basic assumption of our analysis has been the 
use of a simple tunneling form for the intergrain hopping. 
For TFRs, which show a quite weak temperature depen- 



dence of resistivity, this should not be a too serious ap- 
proximation, and simple tunneling captures the essential 
physics, at least as regards piezoresistance effects. For 
other granular metals, or for TFRs at very low temper- 
atures, other processes than intergrain tunneling (grain 
charging effects. Coulomb gap, etc.) play an active role 
in transport. However, as a function of metallic concen- 
tration, we do not expect qualitative deviations from the 
power law behavior of x in the critical region. Instead, 
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